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Abstract. There is an intimate link between supermassive black hole (SMBH) mass (M.) 
and the stellar velocity dispersion (cr.) of the host bulge. This has a fundamental impact 
on our understanding of galaxy and SMBH formation and evolution. However, the scatter, 
slope and zero-point of the relation is a subject of some debate. For any progress to be made 
on this relation, the established values of M. and a, must be robust. Over 50% of current 
M. estimates have been made using the technique of stellar dynamics. However, there is 
serious concern over this method that prompts their re-evaluation. In addition, it is not clear 
how best to define cr,. The aim of the M. - cr. Project is to use STIS long-slit spectroscopy, 
integral field spectroscopy and the latest stellar models, to best estimate the values of M. and 
cr, in as many cases as possible. The project will determine the most appropriate properties 
of the M. - cr, relation itself. 
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1. Introduction 

The discovery of the relationship between su- 
permassive black hole (SMBH) mass (M.) 
and the stellar v elocity dispersion ( c r,) of 
the host bulge (Ferrares e & Merritll |2P00'; 
[Gebhardt et al. 2000) signaled a new and fun- 
damental understanding on the nature of black 
hole and galaxy formation and evolution; the 
M.-cr, relation intimately links the most basic 
characteristic of a SMBH with an underlying 
dynamical property of the surrounding galaxy. 
Consequently, the M, - cr, relation, that is typ- 
ically described by log M, - a+fi log(cr,), has 
received an extraordinary amount of attention. 

The M, - cr, relation is a member of 
a family of scaling relations all linking M, 
to various properties of the host bulge, and 



possibly the entire gravitational mass of the 
host galaxy. For a compre hensive review of 
all the scaling relations see iFerrarese &Fordl 
(l2005h . Despite all these other relations, M. - 
cr, has traditionally been the prefeiTed as, at 
its time of discovery, it showed a remarkably 
small scatter; a smaller scatter than for the 
other scaling relations and a scatter consistent 
with the measurement errors alone. Since then, 
howev er, Marconi & Hunt (2003) and Graharr] 
(I2007h have shown that with enough care, scal- 
ing relation scatters can be comparable espe- 
cially when using near-infrared luminosities. 
Regardless of which scaling relation is "best", 
it is important to fully investigate them all as 
the intrinsic scatters themselves contain pivotal 
information on the precise inter-play between 
SMBHs and their hosts. 
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Reproducing the features of the M, - o", 
relation has become an integral part of SMBH 
and galaxy formation and evolution models. 
The relation also has several other significant 
applications. For example, the M, - cr^ rela- 
tion may be extrapolated to smaller and larger 
values of M,, and be used to estimate M, in 
systems where the black hole sphere of influ- 
ence radius (r/,) cannot be resolved. It is there- 
fore especially important to know the values of 
a and j3, and their associated scatters, to high 
accuracy. To have doubt in either parameter is 
to introduce large uncertainties in extrapolated 
values of M,. This in turn has a fundamental 
impact on our understanding of both black hole 
and galaxy formation and evolution. 



2. The M, - cr^ Relation 

The M, - cr, relation has had many attempts to 
fit its slope, zero-po int and scatter to the whole 
galaxy populat ion ( Ferrarese & Merritt 20001: 
Gebhai-dt et al.' '2000 t iMerritt & Ferraresg 
2001; Tremaine et al. I2002h in addition to 



individual families of gala xies that potentially 
lie off the r elation (e.g. IWang &Lul 1200 lb 
iMathur & Gru pe 2005; Ko mossa & xi1l2007l) . 
So far, these individual fits have produced 
considerably different results (J3 ranges from 
3.8 to 4.9, for example). To demonstrate what 
affect these uncertainties have on the impor- 
tance of the M, - cr, relation, we have used the 
direct M, and cr„ data compiled recently by 
dFerrarese & Fordll2005L Table 2) to constr uct 
Figure [B The fit of lTremaine etalj (|2002|) is 
also plotted as these values are typically used 
by many authors. Taking cr, - 20km s"' (a 
large globular cluster) the two slopes predict 
logM. = 4.IIM0 or log M. = 3.36M0. Taking 
(T, = 470km s"' (a brightest cluster galaxy) 
the two slopes predict logM, - 9.62Mq or 
logM. = 1O.O2M0. 

The aim of the M, - cr^ Project is to pro- 
duce the most reliable estimate of the M, - cr^ 
relation to date. It will couple ground-based in- 
tegral field spectroscopy (IFS) with high spa- 
tial and spectral resolution data, from the STIS 
archive, to estimate M, in a large sample of 
galaxies using the latest three-integral, axisym- 
metric, orbit-based stellar dynamical models. 




Fig. 1. The M, - cr^ relation (where r/, is re- 
solved) using thei erra rese & Ford (2005) {[5 — 
4.86, solid line) and .Tremaine et al.i (i2002h 
OS = 4.02, dashed line) fits. 



The project will also examine the accurate de- 
termination of cr, within the host bulges. In 
addition, the project will identify and address 
several key issues that surround the uncertain- 
ties in determining a and fi. 

2.1. Estimating M, 

Four methods have been used to directly esti- 
mate the values of M, defining the M, - cr, 
relation. Prope r motions have been used in 
the Milky Way dOenzel et alJll997HGhez et alJ 
12005) and H2O masers have be en used to esti - 
mate M. in NGC 4258 (Mivoshi et alJll995l) . 
All other estimates have been made using 
gas and stellar dynamics. Gas dynamics (e.g., 
Ferrarese et al. 1996; Marconi et al. 2001), are 
relatively easy to model as long as the nuclear 
disk is completely dominated by the gravita- 
tiona l potential of the SMBH. Stellar dynamics 
(e.g., van der Marel 1994;' Verolme et alJ2 00?. 
Gebhardt et al. 2003), does not suffer from this 
problem, but is a more complex method that re- 
quires a high sensitivity instrument. It is diffi- 
cult to judge the consistency between estimates 
made by each method as very few galaxies 
have been modeled by more than one method. 

Given a sufficiently sensitive instrument, 
and a sufficiently elaborate model, then stel- 
lar dynamics is the preferred method for de- 
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termining M,. Indeed, to date 57% of direct 
M, estimations have been made using stellar 
dynamics; it is the d ominant method defi ning 
the M. - cr^ relation. ' Valluri et all ((2004') pro- 
vide a detailed description of an orbit-based 
model designed to recover the gravitational po- 
tential of a system using stellar kinematics. 
They carry out detailed analytical testing of 
this model a s applied to a dataset from M32. 
IValluri et al.l ( 12005) apply the same model to 
constrain M. in NGC 205. It is found that a 
large range in M, gives equally acceptable fits. 
The degeneracy in the solution is a function of 
the number of orbits used in the models and the 
number of constraining data points. 

As M32 and NGC 205 are the only datasets 
for which space in stellar dynamical mod- 
els has been explored in detail, it is clear that 
the dominant mass estimates within the current 
M, - o", relation may include un-investigated 
issues that could have a fundamental impact 
on our understanding of the relation itself. The 
M, - cr, project will address these issues in 
two key ways. Firstly, it will use a wide range 
of constraining data points by employing IFS. 
Secondly, it will serially reduce a large family 
of dynamical models, that have a wide range 
of total orbit numbers, on a high performance 
computing cluster with a large node number. 

It has been r e peated by many a uthors 
dPerraresd [ 2002t iMarconi & Hun3 120031: 
IValluri et al.l l2004 that resolving ru (given 
by GM./crl) is a minimum requirement for 
detecting a SMBH. If r/, is un-resolved then 
estimates of M, can be tarnished by uncertain 
contributions from stellar masses. The deriva- 
tion of r/, was first attempted by Peebles ( 1972) 
based on the ad hoc assumption that in a steady 
state the stellar density is exponentially linked 
to the total gravitational potential and that 
o", is constant. The black hole removes stars 
at small radii (via accretion or ejection, for 
example) and disrupts this equilibrium. The 
steady state is restored via an increase in cr» of 
stars brought in from larger radii. 

Although the requirement to resolve r/, has 
been thoroughly argued it must be noted that 
this creates a somewhat degenerate problem; 
you are forced to make an assumption about 
SMBH demographics (i.e., that they follow the 




o» (km s-') 

Fig. 2. The ability of HST to resolve r/,. Above 
the dashed lines r/, is resolved and M, esti- 
mates can be made. 



M, - (Tf relation) in order to deduce whether 
you can detect that SMBH in the first place. 
Furthermore, you are using the very parame- 
ters in the M, - cr, relation in the calculation 
of r/,. Nevertheless, with the assumption that 
SMBHs do follow the M, - cr, relation, we can 
plot ri, within the M,-cr^ plane. In addition, we 
can determine the ability of HST, from which 
most M, estimates are derived, to resolve r/, at 
a specific distance (Fig|2]i. 

Two key points are associated with 
Figure |2l Firstly, it does not include the sen- 
sitivity, and therefore the SMBH detection ef- 
ficiency, of the telescope. The relatively small 
aperture of HST ensures that many STIS orbits 
are required before a sufficiently large contin- 
uum signal-to-noise is reached for stellar dy- 
namical studies. Secondly, it shows that even at 
16 Mpc HST is unable to explore a significant 
area of the M, - cr^ plane (below the dashed 
lines). If there were a population of bona-fide 
under-ma s sive b lack holes, as suggested by 
IVolonteril (|2007|) for example, HST would not 
be able to detect them; they would not be mod- 
eled and the M, - cr, relation would still re- 
main to appear universal. In short, HST has a 
SMBH "blind spot". To unequivocally define 
the M, - cr, relation the entire M, - cr^ plane 
needs to be detectable over a large enough 
distance to include a large population of host 
galaxies. 



4 



Batcheldor: The M, - a. Project 



2.2. Estimating cr^ 

One of the problems associated with the fitting 
of the M, - Q-f relatio n concerns the values 
of cr» used. iFerrarese & Merritt (2000) trans- 
formed cr, to an aperture of r^/S, w here is 
the ho st eifective radius, while Gebhardt et alj 
(120001) use a value weighted by the host lu- 
minosity. The use of different values is un- 
derstandable considering there is, as yet, no 
clearly defined aperture size through which cr» 
should be determined. 

When considering the nature of a LOSVD, 
the second order of which is usually described 
as (Tt, one can see that the size and nature 
of the aperture through which the LOSVD is 
collected has a fundamental bearing on its ob- 
served shape. In addition to t he intrinsic varia- 
tion o f cr* with aperture size dBatcheldor et alj 
l2005h . a velocity field across the aperture will 
rotationally broaden the LOSVD as a function 
of the aperture diameter. This will produce an 
over estimate of the intrinsic cr». In the case 
of an axisymmetric velocity field, the aperture 
shape and orientation will also have a bearing; 
a different value of cr, will be derived from, for 
example, a long-slit aligned along the major or 
minor axis of rotation. IFS can be used to ad- 
dress these concerns. 



3. Summary 

The M, - cr» relation is one of the most valu- 
able correlations in modern galaxy and SMBH 
formation and evolutionary models. There are, 
however, concerns in the ways both cr, and 
many M, estimates have been determined. The 
latest stellar dynamical models have shown 
that great care must be used to fully explore 
space around M,. It is essential that r/, is 
fully resolved, that there are large number of 
high quality constraints, and that a large vari- 
ety of orbit families are employed in the mod- 
els. The most appropriate values of cr, may not 
have been determined as it has yet to be clearly 
defined. STIS long-slit, and integral field spec- 
troscopy can be used to address all these issues 
and will lead to the most accurate form of the 
M, - (7^ relation to date. 
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